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Abstract. Using the persistent current I induced by an Aharonov-Bohm flux in square lattices with random 
potentials, we study the interplay between electronic correlations and disorder upon the ground state (GS) 
of a few polarized electrons (spinless fermions) with Coulomb repulsion. K being the total momentum, 
we show that I oc K in the continuum limit. We use this relation to distinguish between the continuum 
regimes, where the lattice GS behaves as in the continuum limit and I is independent of the interaction 
strength U when K is conserved, and the lattice regimes where I decays as U increases. Changing the 
disorder strength W and U, we obtain many regimes which we study using the map of local currents carried 
by three spinless fermions. The decays of / characterizing three different lattice regimes are described by 
large U perturbative expansions. In one of them, / forms a stripe of current flowing along the axis of 
the diamagnetic Wigner molecule induced by large electronic correlations. This stripe of current persists 
in the continuum limit. The quantum melting of the diamagnetic molecule gives rise to an intermediate 
"supersolid" regime where a paramagnetic correlated pair co-exists with a third particle, before the total 
melting. The concepts of stripe and of supersolid which we use to describe certain regimes exhibited by 
three spinless fermions are reminiscent of the observations and conjectures done in other fields at the 
thermodynamic limit (stripe for high-Tc cuprates, supersolid for Helium quantum solids). 

PACS. 71.10.-w Theories and models of many-electron systems - 73. 21. La Quantum dots - 73.20.Qt 
Electron solids 

1 Introduction because it forms a Fermi glass with Anderson localiza- 

5-h \ tion at high densities or because it forms a pinned Wigner 

'1-1 Motivations crystal at low densities. Since this transition is observed 

£^ ' at low densities, it was suggested [2] that the low density 

O ■ Let us first give some reasons which have motivated this insulator would be a pinned Wigner crystal, which would 

1 study of a few polarized electrons (spinless fermions) give an unexpected metal when it melts above the critical 

> and of their persistent current in two dimensional dis- density. While a 2d insulator is expected, the observation 

ordered lattice models. In the last decade, the conduc- of a 2d metal remains unexplained. If one goes from 

tance of two dimensional electron gases (2DEGs) have the Fermi glass towards the pinned Wigner crystal via a 

^ been measured [T] at low temperatures as a function 2d metal as the density decreases, the existence of a new 

of their densities. Using 2DEGs created in Ga-As het- metallic phase between two insulating phases of different 

erostructures or Si-MOSFETs at densities smaller than nature becomes unavoidable. And indeed, exact numerical 

those used for studying the weak-localization corrections studies of small disordered clusters have shown [3115111)117] 

to the Drude-Boltzmann conductance, a "two dimensional the trace of an intermediate regime consistent with this 

metal-insulator transition" has been observed. When the hypothesis. 

temperature increases, the conductance does not change Subsequent numerical studies of similar systems have 

at a certain critical density, while it increases at smaller shown 8,9 that an intermediate regime persists with- 

densities and decreases at larger densities. The increase out disorder. The possibility of having some hybrid phase 

corresponds to the insulating phase and the decay to made of a quantum solid co-existing with delocalized de- 

the metallic phase. In Ga-As heterostructures or Si- fects, was suggested in refs. [HUH], in analogy with the "su- 

MOSFETs, there is a more or less important disorder due persolid" proposed long ago by Andreev and Lifshitz for 

to charged impurities and rough interfaces. In presence of quantum Helium solids in three dimensions [1UII1 lj . Us- 

disordcr, a 2DEG can be insulating for two reasons: either ing a similar analogy with Helium physics, but assuming 
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macroscopic inhomogeneities of the carrier density due to 
an inhomogeneous substrate potential, a qualitative expla- 
nation of the dependence of the resistance as a function of 
a parallel magnetic field and of the temperature has been 
discussed in Ref. |T2). 

For a simple 2DEG without a random substrate, to 
understand how the ground state (GS) changes when one 
varies the carrier density remains a challenge despite a 
few decades of efforts. The simplest model for a 2DEG is 
a continuum Hamiltonian H c having three parts: one body 
kinetic terms, a two body Coulomb repulsion and a one 
body potential describing the effect of the positive ions. 
The last potential is necessary to have a stable system of 
electrons because of Coulomb repulsions. However, if one 
takes periodic boundary conditions (BCs) and a uniform 
jcllium for the positive ions, the corresponding potential 
becomes a constant term. The scale of the one body quan- 
tum effects is given by the Bohr radius as = fi 2 /me 2 , e 
and m being the electronic charge and mass. The strength 
of the Coulomb repulsion depends on the radius a of a 
circle which encloses on the average one electron. Measur- 
ing the energies in rydbergs (IRy — me 4 /2h 2 ) and the 
lengths in units of a, the continuum 2DEG Hamiltonian 
for N electrons reads 

1 N 2 1 

^ = -^E V ' + 7 E L— —\+ const, (1) 

r * i=l rs l<i<j<N |r>i rjl 

which only depends on a single scaling ratio 

r. = — (2) 
a B 

when N — > oo. We neglect the spin degrees of freedom, 
considering fully polarized electrons (spinless fermions). 
Even in this simpler case, the interplay between Coulomb 
repulsion and the kinetic energy is a complicated issue. 
When many electrons are inside the quantum volume a 2 B , 
r s is small and the 2DEG is a Fermi liquid, with a Fermi 
energy much larger than the Coulomb energy. When r s is 
large, the volume per electron a 2 is large compared to a 2 B , 
and one has an electron solid (Wigner crystal) with weak 
quantum effects. Between those two limits, the nature of 
the ground state remains unclear. 

From Quantum Monte Carlo studies, it is generally be- 
lieved that there is a single first order liquid-solid transi- 
tion. However this result is not free of certain assumptions, 
because of the well known "sign problem" of the Monte 
Carlo methods applied to fermions. One reason to ques- 
tion the existence of a single transition comes from gen- 
eral considerations about Landau theory of phase transi- 
tions and estimates of interface energies. According to Ref. 
|13| . a single first order transition is ruled out because a 
macroscopic phase separation between a liquid and a solid 
is unstable. A way to fix the sign problem of the Monte 
Carlo approaches consists in imposing the GS nodal struc- 
ture. Using a fixed node approach, two nodal structures (a 
Slater determinant of plane waves for the liquid, of local- 
ized orbitals for the solid) have been compared |14j , giving 
rY ~ 37 for the critical density at which the transition oc- 
curs between the two assumed nodal structures. However, 



one cannot exclude the existence of better nodal struc- 
tures giving lower GS energies when r s is neither small 
nor large. And indeed, a third nodal structure has been 
recently studied fTMTT)] in a fixed node Monte Carlo ap- 
proach, using Bloch waves in a variational attractive po- 
tential located at the sites of the Wigner lattice instead 
of free plane waves. This third nodal structure gives a 
lower GS energy for 30 < r s < 80 and N — > oo, support- 
ing the existence of a new hybrid liquid-solid phase. Since 
one does not know the nodal structures, the Monte Carlo 
methods make problem, and exact numerical studies of 
small systems remain useful, despite the presence of large 
finite size effects. They could provide a better understand- 
ing of the low energy physics of a 2DEG, suggesting better 
trial wave functions for fixed node Monte Carlo methods. 
Or they could be the starting point of a finite size scal- 
ing theory ^7] for the interacting 2DEG. This is the road 
which we have explored in a series of works pllM7ll8ll9lll8l 
119) and which we continue here. 

To do an exact study, we take N electrons on a L x L 
lattice for having an Hilbert space of finite size. If A~ and 
L are small enough, the GS can be obtained by Lanczos 
algorithm. Beside the finite size effects due to the finite 
value of N, one has the lattice effects due to the finite 
value of L. If the lattice effects are irrelevant, the lattice 
GS behaves as the continuum GS. Notably, the parameter 
r s remains a scaling parameter. If the lattice effects are 
relevant, r s becomes meaningless, the scaling behavior of 
the continuum limit ceases to exist for giving rise to a new 
lattice physics. In the absence of disorder, the lattice ef- 
fects were studied ^U] as a function of the strength U of 
the Coulomb repulsion. The strength U* above which the 
lattice GS and the continuum GS become different was 
given. The previous studies of a few particle systems [HI 
give the following picture: If one increases U for a 
sufficiently small number N, the Fermi system disappears 
above U F to become a Wigner solid above U w , through 
an intermediate regime for U F < U < U w which could 
be the microscopic trace of an "electron supersolid" , be- 
fore exhibiting strong lattice effects above U* . The Wigner 
solid without lattice effects can be described by a contin- 
uum expansion where the small parameter is 1/ a/tJ. For 
a lattice model, r s = {UL/t){\/\jAnN), where the hop- 
ping term t sets the scale of the lattice bandwidth without 
interaction. Above U* , the continuum expansion breaks 
down and the Wigner solid with lattice effects is described 
by a different expansion in power of another small param- 
eter t/U. If A~ is large, the lattice effects can appear in 
the Fermi regime, a case which we will not consider here. 
In Ref. j the different regimes have been illustrated by 
numerical studies using A~ = 3 spinless fermions in lat- 
tices of size L — 6 up to 18. Three criteria where given to 
get U* and check for A" = 3. The numerical results have 
been reproduced by a continuum expansion of the 

Wigner solid for U < U* and by a lattice i//7-expansion 
for U > U* . Since it is possible to extend these expansions 
to an arbitrary number N of particles, the lattice threshold 
r* below which the lattice effects are relevant was given 
in the thermodynamic limit as a function of s/ob, s being 
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the lattice spacing. The purpose of this work is to use the 
persistent current as a tool for characterizing the effect 
of a weak disorder upon the four regimes which we have 
identified without disorder. The considered system is a 2d 
torus pierced by an Aharonov-Bohm flux. The random po- 
tentials remove the various translational, rotational, and 
inversion symmetries exhibited by the non-disordered lat- 
tice. The persistent currents display new regimes induced 
by the interplay between electronic correlations and ran- 
dom potentials and show more precisely the nature of the 
supersolid GS. In the numerical part of this work, we take 
N = 3 spinlcss fcrmions for lattices of sizes L — 6 and 9. 
When possible, our numerical results are reproduced by 
analytical expansions which can be generalized to more 
particles. 

1.2 Summary of the main results 

This paper is organized as follows: the disordered lattice 
model is introduced in sectional The persistent currents 
(local, total, transverse, longitudinal) yielded by an en- 
closed Aharonov-Bohm flux <P are defined in section [3] A 
relation between the expectation value (K) of the total 
momentum and the current / is derived in section 0] for 
the continuum limit of a lattice of arbitrary dimension 
d. Since this relation is only valid when the one particle 
states of high momenta k are empty, the continuum and 
the lattice regimes can be defined from the study of /. In 
the continuum regime, / is invariant when K is invariant. 
In the lattice regime, / can decay while K remains invari- 
ant. In the continuum regime, the lattice GS behaves as 
in the continuum limit and displays the same universal 
scaling laws. In the lattice regime, the lattice GS becomes 
different and there is no scaling. 

Hereafter, we assume N — 3 spinlcss fcrmions and we 
vary the lattice parameters t (kinetic energy), U (Coulomb 
repulsion), W (disorder) and L (lattice size). In section 
we review the 4 regimes obtained when U/t increases 
without disorder (W/t — 0) for a sufficient value of L: 
continuum Fermi system (U < U F ), continuum super- 
solid (U F < U < U w ), continuum Wi gner molecule 
(U w < U < U*) and lattice Wigner molecule (U* < U). 
The effects of a very weak disorder are described in sec- 
tion One gets 3 lattice regimes as U exceeds U*. First 
the disorder remains negligible and one gets a Ballistic 
Wigner Molecule (BWM regime) . When U reaches a value 
U stripe > U* , a new lattice regime characterized by a 
Coulomb Guided Stripe of Current (CGSC) is found when 
U stripe < U < Ui oc . Above Ui oc , one gets a Localized 
Wigner Molecule (LWM regime). In the CGSC regime, / 
flows along the axis of the Wigner molecule instead of flow- 
ing along the shortest direction enclosing the flux, as when 
U < Ustripe (BWM regime) or when U > U loc (LWM 
regime). The perturbation theories describing these three 
lattice regimes are given, and their range of validity are 
estimated. In section [3 figures where one can see how / 
depends on U, W and t are shown, exhibiting the decay 
of I predicted by the three lattice perturbation theories, 
and the continuum regimes where / is invariant when K 



is invariant. In section |S1 the effect of the disorder upon 
the continuum-lattice crossover is studied. We also show 
how an avoided level crossing induces a jump of I and a 
change of its sign when N = 3, in a continuum regime 
where / is otherwise invariant. The phase diagram of the 
different lattice regimes obtained using N — 3 spinless 
fermions is sketched in the plane (U/t, W/t) fFig llSJI . In 
section|21 we study the effect of random potentials in con- 
tinuum regimes where I remains invariant. We give a de- 
tailed study of the case W = t = 1, where the motion 
is diffusive without interaction for the considered values 
of L. The study of the map of local currents shows that 
the stripe of current observed in the lattice CGSC-regime 
persists in the continuum limit. The persistence of Id mo- 
tions yielded by the Coulomb repulsion in a 2d disordered 
lattice allows us to obtain the parity of the number of 
particles contributing to I from the sign of /. This leads 
us to suggest that a A^ = 3 diamagnetic Wigner molecule 
melts in a disordered lattice through an intermediate "su- 
persolid" regime where a. N — 2 paramagnetic Wigner 
molecule co-exists with a third nearly localized particle, 
before becoming a Fermi glass. In sectional the three first 
harmonics of the function I(^>) are given as a function 
of r s , showing in greater details how a 2d current I of 
random sign becomes a Id current of given sign when the 
electronic correlations increase. In sectional we underline 
the main results obtained in studying three particles: the 
existence of a regime of stripe and of a supersolid regime 
where the zero point motion of the electron solid becomes 
of the order of the system size without disorder, or where 
a delocalized pair co-exists with a localized particle with 
disorder. We conclude by discussing the possible relevance 
of our results when A~ becomes larger. 

2 Lattice Hamiltonian with disorder 

The lattice Hamiltonian Hi describing A" polarized elec- 
trons free to move on a L x L disordered square lattice 
with periodic boundary conditions (BCs) reads: 

Hi = Hkin + Hint + Hdis- (3) 

We write these three terms using the creation (annihila- 
tion) operators in real space and in momentum space. 

2.1 Site basis 

Using the operators cj (cj) which create (annihilate) a 
polarized electron (spinless fermion) at the lattice site j, 
the kinetic term reads: 

H kin = t I AN - c h' ' ( 4 ) 
\ UJ'J / 

the Coulomb term reads: 

Hint - ]X7p ^ 
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The third term describes a random substrate: 



Ha 



(6) 



W gives the strength of the disorder and the variables 
€j are taken at random between [—1/2,1/2]. means 

that the sum is restricted to nearest neighbors, nj = cjcj. 
For a square of area D 2 , the lattice size is given by L — 
D/s, where s is the lattice spacing. Assuming periodic 
BCs, we define the distance d$y between the sites j and j' 
in unit of s, as 



d ii' = 



L 



sin 



L 



2 fek 

L 



(7) 



This metric coincides at short distance with the natural 
2d metric and avoids |19j to have singular repulsions with 
cusps when djy k, L/2 and s — > 0. The corresponding 
Coulomb repulsion is essentially equivalent to the one ob- 
tained from Ewald's summation, where one assumes the 
infinite periodic repetition of the same L x L square in- 
stead of a 2d torus. 



2.2 Plane wave basis 

The Hamiltonian J3J can also be written using the oper- 
ators d\, (dk) creating (annihilating) a polarized electron 
in a plane wave state of momentum k: 



T.T, 



(8) 



K, 



Ki 



H 



int 



H 



dis 



Fig. 1. Structures of the interaction and disorder Hamiltonian 
matrices written in the basis of the PWSDs. The kinetic en- 
ergy gives a diagonal matrix. We have sketched the structures 
for 4 subspaces of total momenta Ki, . . . , K4. The matrix ele- 
ments are zero in the white parts. For Hi nt , the darker regions 
corresponds to hopping terms between PWSDs built out from 
one body states of small momenta k only. In the shaded parts 
of Hint and Hw, the matrices are sparse. 



all the components of k are small such that sin k x rj k x 
and sin k y « k yi are indicated by a darker grey. Moreover, 
only PWSDs of same K having N —2 momenta k in com- 
mon out of N can be directly coupled by the two-body 
interaction. When N > 3, this means that each block of 
given K is sparse. In contrast, the random matrix due to 
H\y has zero matrix elements inside the diagonal blocks 
of momenta K, but gives non zero off-diagonal terms cou- 
pling those diagonal blocks. These coupling off-diagonal 
blocks are even sparser, having non zero terms between 
PWSDs differing by a single k only, out of N. 



Hh 



ANt - 



2t (cos k x 

k 



cos ky) dj^cik (9) 



flint = U J2 V (^) d 



L+ q 4'-<A'dk 



k,k' 



where 



and 



ffdia=W^e(q)4 +q d k 

k,q 

1 cosc U 



3j0 



— ^ejexpHq-j). 
j 



(10) 



(11) 



(12) 



(13) 



In Fig. ^ the structure of the Hamiltonian matrix for 
N spinless fermions is sketched in the basis of the Plane 
Wave Slater Determinants (PWSDs) <4 . . . <4 N |0>. The 
matrix corresponding to H^m is diagonal. If we order the 
PWSDs by series of same total momentum K = J^iLi 
the matrix corresponding to Hi nt is block-diagonal, each 
block being characterized by the same momentum K. In- 
side each block, we order the PWSDs by increasing kinetic 
energy. The matrix elements between PWSDs for which 



2.3 Lattice factor r s 

The continuum thermodynamic limit is obtained when the 
finite size effects due to the finite value of N and the lattice 
effects due to the finite value of L are negligible. This is 
clearly out of reach of the Lanczos algorithm. However, for 
a finite value of N, r s remains a scaling parameter when 
the lattice effects are negligible. In a medium of dielectric 
constant e, the strength U of the Coulomb repulsion reads: 



U = 



zS 



while 



t = 



(14) 



(15) 



2m* s 

is the hopping element between nearest neighbor sites and 
sets the bandwidth of carriers of effective mass m* . Taking 
to = to* and e = 1, one has in lattice units lRy = U 2 /At, 
as = 2st/U and 



UL 



a B 2WttN 

for the rydberg, the Bohr radius and the factor r s 
tively. 



(16) 



respec- 
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Instead of changing the number N of particles, one can 
vary r s keeping N fixed and varying U, t and L, and hence 
the parameter 



n 



UL 



= rsVlnN. 



(17) 



This amounts to change r s by changing as instead of a. 
The lattice effects are important ^5] either at large values 
of U (hence large 77) for a given TV, or at large values of 
N (hence small r s ) for a given U. Hereafter, we use n or 
U instead of r, in order to avoid confusion. 



The latter equation comes from the conservation of the 
current. The total transverse current It(&) is given by 
a similar formula using the 

The change AEq (<?) of the GS energy E induced by 
<Z>: 

AE Q ($) = £ (0) - E ($) (21) 
is a quantity related to I which is simpler to study. 



The GS total longitudinal current 



t 27Tin 



(22) 



3 Persistent current / 

An Ahronov Bohm flux <P is enclosed around the longitu- 
dinal x-direction as sketched in Fig. [21 



can be expressed in terms of the operators d\, (d*.)- Here- 
after, the k x and k y components are quantized in multiples 
of 2ir/L, and an explicit phase factor exp(2ir<P/(<PoL)) is 
added to the ^-components. This gives 



(23) 







4 



















I 



Fig. 2. 2d Torus with N electrons enclosing an Aharonov- 
Bohm flux $ oc B 



By a gauge transformation, <P can be included in the 
longitudinal BC (antiperiodic BC corresponding to <P = 
<Pq/2 in our convention), while the BC remains always 
periodic in the transverse y-direction. This flux creates 
persistent currents for a GS of wave-function |<?o) an d of 
energy E ($): 

— At a site j = (j x ,jy), the local longitudinal current 
j- (<P) and the local transverse current jj {&) read 



jU<P)=2%(% 



#0) , (18) 
(19) 



in units of \Jtj2mj s. 

The total longitudinal current Ii(<P) enclosing the flux 
is obtained by summing the local longitudinal currents: 



dE Q (<P) 
d<P 



.),,= 



j 



(20) 



for the total longitudinal current and 



(24) 



for the total transverse current. 

Expressed in terms of its projections •Z'ki k N onto the 



PWSDs n„=i dl n |0), the GS reads: 



l*b) 



E ^i,-.kN II 0> (25) 

ki....,kisj n—1 



which gives 



.,k> 



i=l 
N 



it($) = I E ^...k.i'E 8111 ^- 



(26) 



ki, 



( = 1 



4 Relation between 7 and K in the 
continuum limit and lattice threshold U* 

As far as the lattice GS occupies small momenta k, one can 
approximate the dispersion relation of the lattice model by 
the parabolic dispersion relation of the continuum limit: 

2 — 2 cos ki w kf and sin ki fa ki. (27) 

In this approximation, the expressions H26JI become 



N 



L 



E '".<, kx 1 E A> 



hi'J'^l E l*ki,...,kN| a E**: 



1=1 

iV 



2tt<2> 



(28) 
(29) 
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which give 



2nN $ 
L <? 



h = T (Ky) 



(30) 
(31) 



where (K XtV ) = E kl ,...,k N (E 4 =i h<x, V) 

This proves an important theorem which was known 
in one dimension and numerically checked 20,21,22] in 
two dimensions: the GS persistent current I is indepen- 
dent of the interaction strength U in the continuum limit 
of a non-disordered lattice model of arbitrary dimensions. 
Moreover, this provides one criterion (criterion 2 of Ref. 
|19|1 for obtaining the lattice threshold U* . We have shown 
that / is proportional to K when the momenta k occupied 
by the GS are small, allowing to use the approximation 
(|27() . The invariance of J is a direct consequence of the 
invariance of K. When U = W — 0, the approximations 
(|27|l require small filling factors N/L 2 . If one turns on U 
for W — and a small filling factor, J a K remains in- 
dependent of U when U is smaller than the threshold U* . 
When U > U*, the GS of a non-disordered lattice begins 
to occupy, inside the same subspace of total momentum 
K, large one particle momenta k, and the relation be- 
tween / and K ceases to be valid. Without disorder, K 
is always conserved, and I remains independent of U as 
far as I cx K. The breakdown of the invariance of / is a 
consequence of a localization-delocalization crossover in- 
duced by U in momentum space. For a fixed value of t, this 
crossover can be induced by U inside a subspace of given 
total momentum K, by W among subspaces of different 
K, or by the combined role of U and W. The regime where 
K is conserved defines the ballistic many body regime. 



5 Overview of the non disordered case for 

N = 3 

Before considering the case where random potentials are 
included, it is useful to summarize what we know when 
there is no disorder. Let us review a few published 9 
IT§] and unpublished results based on a study of the case 
N = 3. 



5.1 Continuum Wigner molecule above U w 

(r s > rf re 50) 

On a continuum D x D square with periodic BCs, the 
Coulomb repulsion with the taken 2d metric reads 



V{v) 



£ A /sin 2 ^ + sin 2 ^f 



(32) 



The GS center of mass is delocalized. When the den- 
sity N/D 2 is small enough, the spacings d between the 
particles are large and one obtains a Wigner molecule, 



which is in the continuum limit if the fluctuations of d 
are also larger than s. For a given center of mass, to put 
the particle on the sites ri = (0,0), T2 = (D/3, D/3) and 
r3 = (— -D/3, —D/3) yields for the Coulomb energy a min- 
imum value: 



E, 



Coul 



V6e 2 ir 
D 



(33) 



For a large D, one can expand the pair-potential around 
the equilibrium inter-particle spacing r = {D/3, D/3) up 
to the second order to get harmonic oscillations around 
the equilibrium positions. The Hamiltonian H c pa Ecoui + 



Hharm, where the harmonic part reads: 



H, 



1 ^ 



harm 



2m 



i=l 



XMX. 



(34) 



The vector X = (x\, y\, X2, 2/2> %3> Vs) describes the mo- 
tions of the molecule around equilibrium. The 6x6 matrix 
M is given by: 



M 



where 



/ 2A 


2B 


-A 


-B 


-A 


-B\ 


2B 


2A 


-B 


-A 


-B 


-A 


-A 


-B 


2A 


2B 


-A 


-B 


-B 


-A 


2B 


2A 


-B 


-A 


-A 


-B 


-A 


-B 


2A 


2B 


\-B 


-A 


-B 


-A 


2B 


2A ) 



(35) 



A 



7y/6 



9 

en 



72 (36) 

and B = 3 A/7. Diagonalizing M, one obtains the normal 
modes, while the eigenvalues of M give their frequencies. 

One obtains two eigenvectors Xi and \2 of eigenvalue 
0, which corresponds to the translation of the center of 
mass of the molecule, two other eigenvectors X3 an d Xi 
of eigenvalue 10B, which corresponds to the longitudinal 
mode, while the two last eigenvectors X5 an d X6 °f eigen- 
value AB correspond to the transverse mode. The expres- 
sions of eigenvectors Xi are given in Ref. . 

Using these normal coordinates, the Hamiltonian (|34|l 
becomes a decoupled sum of four harmonic oscillators: 



& 



Hharm „ 7 , o > 

2m ^ d X i 



10B( X l+ X 2 )+4B(xt+xl), (37) 



Without disorder, the interaction can only couple 
states of same total momentum K. Inside a subspace of a 
total momentum K, the wave functions can be factorized 
as 

. . . )X8 ) = ^(xi,X2)Cr"" 6 (X3, • . • ,Xe). (38) 

The kinetic energy associated to the motion of the center 
mass reads: 



£cm(K) 



2iVm 



-K 



(39) 



The wave-function ^"3-°>— >"e— (x3, . . . ,Xe) correspond- 
ing to the relative motions can be factorized as: 



U3— 0....,ri6— 



^0L(X3)^0L(X4)^0T(X5)^0T(X6) (40) 



Z. A. Nemeth and J.-L. Pichard: Persistent currents 



7 



for the ground state (713 = 0, . . . , ng = 0) of momentum 
K. L,T refer to the transverse and longitudinal modes 
and ipo is the GS of an harmonic oscillator: 



(p (x) 



1 



of characteristic length: 



1/4 



(41) 



(42) 



For the expanded pair potentials, the GS energy Eq(K) 
inside a subspace of momentum K becomes 



E {K) = E, 



Coul 



E cm (K) + h(uj T + lu l ) 



where 



20B 



In lattice units, one gets: 



-E-o(K) — Ecoui + E cm (K) + E re ifi, 
where E re i(n 3 = 0, . . . , n 6 = 0) = E relfl . 
t 



S cm (K) 



E, 



N 

V6ttLT 



(43) 
(44) 

(45) 
(46) 



rel,0 



E'o(K) depends on K only through E cm (K), while 
Ecoui and the transverse and longitudinal modes remain 
independent of K. This is still the case if one goes beyond 
the harmonic approximation for the pair potentials, in- 
cluding anharmonic corrections. The only approximation 
which we have done is to assume that the relative motions 
are small compared to D. Thus the GS relative wave func- 
tions !f" e 3 z =0 " " : " 6= °(X3, ■ • • ,Xe) and E re i(n 3 = 0, . . . ,n 6 = 
0) are approximated by the solutions of an unbounded 
2d system, while the wave functions of the center of mass 
^crnixi, Xa) and E cm (K) correspond to a bounded 2d sys- 
tem. The assumed 2d harmonic oscillations make sense 
only if the Wigner molecule is rigid enough for making 
negligible the boundary effects upon the relative motions. 



5.2 Supersolid molecule below U w 

(rf < r s < rf ) 

To determine the values of r s where the GS energy is 
described by Eq (|431) . we plot in Fig. |31 the GS ener- 
gies E (K) - E cm (K) - E C oui as a function of r s , for 
a size L — 18 and different subspaces of momentum K. 
Ecoui = V6nU/L when L/3 is integer. 

The solid line (0.139-^/7-7) gives the behavior implied 
by the continuum harmonic expansion fEa I45|) with un- 
bounded relative motions. One can see that this theory is 
accurate for r. < r s < r* Above r. ~ 50, the relative 



10 

3 
o 
O 
LU 

E 

° -I 
LU 1 

o 
LU 



0.1 




K=(0,0) 
K=2ji/L(1,1) 
K=2ji/L(2,1) 
K=2ji/L(2,2) 
K=2ji/L(3,0) 
K=2ji/L(3,3) 



0.1 



10 100 1000 10000100000 



To 



Fig. 3. W = 0, L = 18 and N = 3: The GS-energies E (K) - 
i?cm(K) — Ecoui for different total momenta K as a function 
of r s . The arrow shows rY ~ 50 above which the behavior is 
given by Ea 1451 fsolid line 0.139-^/?!) up to the lattice threshold 



motions become localized and coincide to unbounded har- 
monic vibrations. Below r^, the behavior is more compli- 
cated. On one side, the Eo(K) are roughly described above 
rf « 10 by the expansion (|45|l valid for a continuum solid 
molecule. On the other side, the relative motions depend 
on K, a dependence which cannot exist unless the relative 
motions are extended over a scale of the order of the sys- 
tem size D — Ls. Large zero-point motions are expected 
in a supersolid. For TV = 3, we have shown that there is 
regime where there is a floppy Wigner molecule, though 
the relative motions depends on the quantization of K and 
hence on the BCs. This intermediate "supersolid" regime 
was noted in Ref. [5] using another 2d metric for L = 6, a 
limit where the lattice effects play a role. Our results for a 
Wigner molecule describable by a continuum theory when 
L = 18 gives further support to the existence of an inter- 



mediate regime for 



< r. < r 



w 



which is not a lattice 



effect. Hereafter, we use the word "supersolid" to refer to 
this regime. 



5.3 Lattice regime above U*(L) 

(ri>r*(L)) 

The GS energy of the lattice model (Hamiltonian |2Jl is 
also not described by the continuum expansion (|45|l above 
the lattice threshold r*(L). Three equivalent criteria were 
introduced in Ref. [T§: for giving r\ (L) . One was based on 
the breakdown of the invariance of I discussed in section 
0] Let us introduce a fourth criterion which gives a similar 
answer: the lattice spacing s becomes relevant when the 
size of the harmonic oscillations becomes of the order 
of the lattice spacing s. The longitudinal modes giving the 
smallest scale, the criterion reads 



-i 2-, ,2 



1/4 



S, 



(48) 
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giving r\ 



(L) « (2V6)/(5^)L 4 



U*{L) 2a/6 



t 



5 71 



L 



(49) 



When U > U*(L) (n > r*(L)), the lattice model exhibits 
a different behavior than its continuum limit and r s ceases 
to be a scaling parameter. This difference is shown in Fig. 
01 where one can see the universal scaling regime below 
r ; * (L) and the non universal lattice behaviors above r\ (L). 

To show scaling when the GSs of different total 
momenta K correspond to a rigid molecule with K- 
independent harmonic oscillations and a K-dependent 
translation of the center of mass, we use the dimensionless 
quantum correction i*b(K, rj) to the classical electrostatic 
energy: 



F (K,n) = 



E (K,n) - E cm (K) - E, 



Coul 



£ (K = 0,Z7 = 0) 



(50) 



Fo(K, n) is shown in Fig. fj] as a function of r/ = UL/t for 
different momenta K and N = 3. One can see that the 
data obtained for L = 9 and 18 are on the same universal 
curve below r*(L), the curve depending on K below a 
value 7j corresponding to w 50. The universal K- 
independent scaling curve 0.2326y^7 is obtained assuming 
that E (K = 0, U = 0) R* 8n 2 t/L 2 as in the continuum 
limit: 



V5 + \/2 /V67r 3 



8tt 2 V 3 
yW2 /18\ 1/4 

V 7T y 

0.5764^ = 0.2326Vn. 



(51) 



Fb saturates above r ; *(L) to a value 4Nt/(8ir 2 t/L 2 ) which 
is independent of U. 



100 
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10000100000 



Fig. 4. IT = 0, JV = 3: The dimensionless GS-energies Fo for 
different total momenta K as a function of n = UL/t for L — 9 
and L = 18. The solid line 0.2326^/n comes from Eq H51H . 



6 Effects of a weak disorder in the lattice 
regime of an electron solid 
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Fig. 5. Sketches of the hopping processes characterizing the 
three lattice regimes of a very weakly disordered lattice when 
U > U* for L = 6 and N = 3. From left to right: Ballistic 
Wigner molecule (BWM) when U* < U < U 3 tri P e', Coulomb 
Guided Stripe of Current (CGSC) when U s tri P e < U < Ui oc 
and one particle motion out of a Localized Wigner Molecule 
(LWM) above Ui oc . 



Without interaction (U — 0) one has three regimes 
(ballistic, diffusive and localized) as W increases. Assum- 
ing Drude formula and Born approximation, the dimen- 
sionless conductance g reads |2*3"] : 



96tt- 



N 
L 2 



(52) 



For a Fermi momentum Uf 
mean free path I becomes 



V4ttN/(Ls) 2 , the elastic 



I 



960FN fj_ 
W 



(53) 



If U = 0, = 3, L w 6 — 9, this gives a ballistic motion 
when W/t ~ 0.01 — 0.1, a diffusive motion when W/t w 1, 
and strong Anderson localization for W/t » 10 — 20. 

Let us first study the effect of a weak disorder in the 
simple lattice regime where the three particles form an 
electron solid (Wigner molecule). Since r/ becomes an ir- 
relevant parameter above r ; *, we simply use the param- 
eter U, U* being the interaction strength corresponding 
to ri(L) for the used values of L and t. For U > U* , 
the low energy physics of the Af-body problem can be de- 
scribed by an effective one particle problem. For A^ = 3 
and L/3 integer, the Wigner molecule is oriented along 
one of the two diagonals of the square lattice. The average 
interparticle spacing (d) = (±L/3,±L/3) and the fluctu- 
ations of d are smaller than s above U* . Taking t = 1, a 
fixed value of L and a weak value of W, one gets three 
lattice regimes when one increases U above U* , charac- 
terized by three different lattice perturbation expansions. 
First, there is a Ballistic Wigner Molecule (BWM) on a 
scale L smaller than the elastic mean free path l cm char- 
acterizing the motion of the center of mass. On this small 
scale, the disorder can be neglected, the motion of the 
center of mass remains ballistic and can be described by a 
non-random single particle lattice model with an effective 
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nearest neighbor hopping term t e ft oc t /U . Above and for N = 3 
an interaction Jj r stripe > U* , the effects of disorder be- 
come relevant and must be included in the effective single 
particle model. This gives a new regime where the inter- 
play between disorder and electronic correlations leads to 
a Coulomb Guided Stripe of Current (CGSC), since the 
current / flows along the axis of the Wigner molecule. 
Eventually, the effective model breaks down when U ex- 



l eff 



= 6- 



( 7V6 Uir 
\ 36 



(59) 



ceeds a large threshold Ui oc > U s tri pe > U* , and we have 
the perturbative regime analyzed in refs. [21.25 for a lo- 
calized Wigner molecule (LWM). While the motion in- 
duced by <P is along the shortest paths enclosing <P in the 
BWM and in the LWM regimes, the motion in the in- 
termediate CGSC regime uses a longer path of smaller 
electrostatic cost. 

In Fig. [SJ we have sketched the different hopping pro- 
cesses characterizing these three regimes when L = 6. 
The left figure corresponds to a Ballistic Wigner molecule 
(BWM). The middle figure describes the Coulomb Guided 
Stripe of Current (CGSC) induced by the interplay be- 
tween disorder and electronic correlations. The right fig- 
ure shows a Localized Wigner Molecule (LWM) carrying 
exponentially small one particle currents. 



6.1 Lattice regime for a Ballistic Wigner Molecule 

When L is smaller than the elastic mean path l cm of the 
center of mass of the Wigner molecule, K does not vary 
and one can consider a subspace of total momentum K. 
Assuming L/3 integer and t/U — ► 0, the GS reads 



Including the flux in the longitudinal hopping terms 
t e ffexp (i2nN(p/<PoL), one gets for the change of the GS 
energy 



/ 2ttN<P 
AE (#,K) = 2t eff cos(K x + - cosK x 

(60) 

In Fig. as in the following ones, the maps of local 
currents characterizing a disordered sample and the cor- 
responding site occupation numbers n? given by 



? = (^(t/)|cj Cj |<W)) 



(61) 



l*o(K)) 



A^2e tKi Cj\0), 



(54) 



are shown using the same convention: the arrows between 
neighboring sites give the local currents, the largest lo- 
cal current found in the sample is shown by an arrow of 
length s, while n? = 1 is shown by a disk of diameter s. A 
small value <P — 0.05^o is used for driving the persistent 
currents. 

In Fig. El the map of local currents are given for a 
weakly disordered sample (W = 0.01 and t = 1) with 
an interaction strength U = 300 above U* and the corre- 
sponding site occupation numbers n?. The map of currents 
and the density exhibit an almost perfect translational in- 
variance in this ballistic regime and gives an illustration 
of the BWM motion characterized by t e f / ■ 



where 



°'i C j + (L/3,L/3) C j + (2L/3,2L/3) 



(55) 



and A is a normalization constant. When t = 0, we have 
2L 2 /3 configurations Cj |0) of identical Coulomb energy 
Ecoui- Without disorder, this large degeneracy is removed 
by a hopping term t e ff oc t 3 /U 2 coupling two nearest 

neighbor configurations Cj |0) and Cy |0), as indicated in 
Fig-GDleft. At this order, one gets an effective one particle 
L/3 x L lattice model described by 



H eff (W = 0) = Wt + Ecoui -teffJ2 C l C 'y- ( 56 ) 

cy) 



where the effective hopping term reads 



{P} ri 7 ( E coui - e Pi 



(57) 



P labeling a series of intermediate configurations coupling 
j and j' of Coulomb energy Ep y . This gives 



l eff 



JJN- 



-L 



3JV-3 



(58) 



Fig. 6. Ballistic translation of a Wigner Molecule in the lattice 
regime (BWM regime), obtained in a disordered sample with 
N = 3, L = 9, W = 0.01, t = 1 and U = 300. Left figure: 
map of currents for $ — 0.05^o- Right figure: corresponding 
site occupation numbers n?. 



6.2 Lattice regime for a Coulomb Guided Stripe of 
Current 

As one continues to increase U, the effect of disorder 
cannot be ignored and the Wigner molecule finishes to 
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be pinned in the GS configuration Cj |0) of minimum 
Coulomb energy for which e } is also minimum. As one 
can see in Fig. [7| (right), a complete localization of the 
center of mass of the rigid molecule is achieved for a dis- 
ordered sample with W = t = 1 and U = 1000. The GS 
probability to occupy a site is almost zero outside three 
sites where C\ |0) has a minimum energy. It is likely that 
this complete localization is obtained after a small diffu- 
sive regime. In this pinned regime, two competing motions 
are now relevant for enclosing the flux, which are sketched 
in Fig. (middle and right). The CGSC motion (Fig. 
middle) is larger when U stripe. < U < Ui oc while the LWM 
motion (Fig. [S] right) dominates the longitudinal motion 
when Uioc < U. Let us first study the CGSC motion. The 
effect of disorder can be included by adding random site 
potentials Wuj- } to the effective one particle model (|56|) : 



Heff(W) = H eff (W = 0) + WJ2 "jC/Cj, 



+ ej+(i/3,i/3) + e.j-(L/3,L/3) 



(62) 



where 

Assuming the Hamiltonian I|ri2[l. the GS flux depen- 
dence is obtained at the order 2L/3 of an expansion 
in powers of t e ff/W, corresponding to the (L/3, L/3)- 
translation shown in Fig. middle. One gets for the GS 
energy change 



AE (<2>) cx 



t 2 e y f 3 cos (2n$/$ ) 

, — 2L/3-1 

Vnw 

t 2L L 4L cos (2n$/$ ) 



which gives for the two sizes L which we have studied 



(63) 
(64) 



AE (L = 6, <2>) cx 
AE (L = 9, <P) cx 



W 3 U 8 

t 18 



w 5 w 



cos(27T<P/<Z>o) (65) 
r cos {2n$/<P ) (66) 



respectively. 

The current / forms a stripe which follows the axis of 
the pinned Wigner molecule for enclosing the flux. The 
path is longer than the shortest path, but it uses inter- 
mediate configurations of minimum Coulomb energy. An 
illustration of such a stripe is given in Fig. [7| (left), ob- 
tained using a disordered sample where W = t = 1 and 
U = 1000: the pinned Wigner molecule carries in the 
CGSC regime a current of components /; = It, in con- 
trast to the BWM and LWM regimes where It is negligible 
compared to //. 



6.3 Lattice regime for a Localized Wigner Molecule 

When U becomes very large, the previous CGSC remains 
and gives the transverse current I t , but the main contribu- 
tion to /; is no longer given by the effective Hamiltonian 
l|62|). This is because one particle hops, using L — 1 in- 
termediate states of high Coulomb energy become more 
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Fig. 7. Coulomb Guided Stripe of Current (CGSC) in the 
lattice regime, obtained in a disordered sample with iV = 3, 
L — 9, W = t — 1 and U = 1000. Left figure: map of currents 
for <P = 0.05#o. Right figure: corresponding site occupation 
numbers n? . 



advantageous than the correlated effective hopping using 
2L — 1 states of smaller Coulomb energy. For N — 3, this 
corresponds to one particle hops along the shortest path 
inclosing the flux (see Fig. fright). To numerically get this 
regime for a very weak disorder and L = 9 requires a very 
large value of U/t, where the numerical results become 
inaccurate. Nevertheless, an illustration of this regime us- 
ing a sample with a large disorder (W/t = 20) is given for 
U/t = 1000 in Fig. |H1 with the corresponding GS occupa- 
tion numbers. 



Fig. 8. One particle motions out of a Localized Wigner 
Molecule (LWM) in the lattice regime: Disordered sample with 
N = 3, L = 9, t = 1, W = 20 and U = 1000. Left figure: map 
of currents for $ = 0.05$o- Right figure: corresponding site 
occupation numbers n?. 



The GS energy change AE$ (<P) reads: 

jLr3L-3 

AE ($) cx cos (2ir$/$ ) (67) 

which gives for the two sizes L which we have studied 

AE (L = 6, #) cx lj cos (2^/^) (68) 
t 9 

AE Q {L = 9, <2>) cx cos (2ir<P/<P ) (69) 
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respectively. 



6.4 Crossovers between the different lattice regimes 

To obtain the threshold values where a crossover between 
different lattice regimes occurs, we compare the GS energy 
changes AEq (<£>) of the different lattice regimes (Eqs. I6UI 
ElandEH). 

Increasing W/t for a ratio U/t > U*/t, the BWM- 
CGSC crossover takes place when W reaches a threshold 



value W, 



stripe 



eff 



This gives: 



Wst^ipe k ft j v 



() 



(70) 



The CGSC-LWM crossover takes place when W/t reaches 
a second threshold 



O 

O 
O 
II 

o 




10000 



Fig. 10. L — 9 and t = 1: Amplitude \h($ = O.O50 O )| of the 
total longitudinal current as a function of U for a disordered 
sample with increasing values of W . The thick solid lines give 
the expected lattice decays. 
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oc 
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(71) 7 Numerical results for iV = 3 particles. 



If one increases U/t above U* /t for a given ratio W/t, 
the BWM-CGSC crossover occurs when U/t becomes of 
the order of 

1/2 



u 



stripe 

1 



(72) 



followed by the CGSC-LWM crossover when U/t reaches 



Ulo 



2L-3 

- L + 3 3L + 9 
L L + 3 . 



(73) 



The study of different samples confirms the dependences 
of the effective hopping terms predicted for the three lat- 
tice regimes as a function of the lattice parameters U, W 
and t in the expected ranges of parameters. We show nu- 
merical results for the total currents I\ and It, and their 
ratios It/h obtained using disordered samples, fixing two 
parameters (U, W or t) and varying the third, when a 
small flux <P = 0.05<?o is enclosed. We give the typical be- 
havior obtained using a single sample, without ensemble 
average. The sample to sample fluctuations are negligible 
for weak disorders, and can be more important when W 
becomes large. 




Fig. 9. L = 6 and t = 1: Amplitude |ii(<P = 0.05<? )| of the 
total longitudinal current as a function of U for a disordered 
sample with increasing values of W. The thick solid lines give 
the expected lattice decays. 



7.1 Persistent currents as a function of U 

In Figs. Inland mil the amplitude of the total longitudinal 
current \Ii\ is shown as a function of U for t = 1 and 
various values of W, for L = 6 and 9 respectively. The 
power laws with the predicted exponents characterize the 
decay of the currents (U~ 2 , [7~ 4L / 3 and U X ~~ L ) in the three 
expected lattice regimes. 

In Figs. 1111 and !12l the ratio \h/Ii\ is given as a func- 
tion of U . Since the stripe of current characterizing the 
CGSC-regime yields a ratio \h/h\ = 1, one can see on 
those figures the values of U where the CGSC-regime takes 
place. For small values of U which are below the lattice 
threshold U* , one can see in Figs. l9l and mil that |/;| and 
| it | are essentially independent of U, excepted around a 
value Uic where an avoided level crossing occurs. The value 
of Uic depends on W and t and exhibits small sample to 
sample fluctuations. This avoided crossing yields a change 
of the sign of Ii, a sharp drop of |/;| and a singularity in 
\It/h\ which can be seen in Figs. 1111 and!12l The contin- 
uum behaviors taking place when U < U* , where and 
It are essentially independent of U instead of exhibiting 
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Fig. 11. L — 6 and t = 1: Ratio I t /h for <P 
function of U for increasing values of W. 
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Fig. 12. L = 9 and £ = 1: Ratio 7 t //j for $ = 0.05<£ as a 
function of U and increasing values of W. 



the lattice decays while \h/h\ 

m 



1, are studied in section 



7.2 Persistent currents as a function of W 

Figs. El and El show for L = 6 and 9 respectively how 
the amplitudes |//| (indicated by lines with symbols) and 
\It\ (indicated by the same lines without symbol) vary as 
a function of W for t — 1 and various values of U . When 
U is large (U = 100 — 1000), one can see the expected 
lattice behaviors: the CGSC-decay h oc ry-^L/s-i) for 
W s tripe < W < Wioc followed above W ioc by the IWM- 
regime where // saturates to a small value independent of 
W. When U is smaller (U = 1 — 10), Jj is first almost 
independent of W, before decaying as IT -1 when 1 < 
W < 10 for L = 6 and 9. 

The transverse current It begins to increase as a func- 
tion of W from the zero value characterizing the limit 
W = 0. This increase of It takes place for the values of W 
where // is roughly independent of W . This increase goes 
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Fig. 13. L — 6 and t = 1: Amplitudes \Ii\ (indicated by a line 
with a symbol) and \It\ (indicated by a line without symbol) 
of the total longitudinal and transverse currents as a function 
of W. The data are obtained using a disordered sample with 
$ = 0.05^0- The symbols and lines corresponding to different 
values of U are defined in the figure. The thick solid lines give 
the W~ 3 , IT -1 , W 1 and W 2 decays respectively. 
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Fig. 14. L = 9 and t = 1: Amplitudes /; (indicated by a line 
with a symbol) and J t (indicated by a line without symbol) 
as a function of W. The studied values of U axe represented 
as in FigHSl The thick solid lines give the IT" 5 , W' 1 and W 2 
decays respectively. 



as W 2 for large interactions (U = 100—1000) and as W for 
small interactions (U = 1 — 10), as underlined in Figll3l 
These two different powers can be explained by a pertur- 
bative expansion in powers of W. When W = 0, GSs of 
different total momenta K are degenerate when U is small 
and It oc W. As U increases, there is a level crossing for 
W = above which there is a non-degenerate GS of mo- 
mentum K = and It oc W 2 when U is large. The increase 
of it is followed by a regime where It ~ h when L = 6 or 
exhibits the same decays (W 1 or W~^ 2L ^~ 1 )) than h 
when L — 9. For large U and above Wi oc , It continues to 
decay as in the lattice CGSC-regime (I t oc T T~( 2L / 3 - 1 )) 
while Ii saturates to a IT-independent value. 
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7.3 Persistent currents as a function of t 



the lattice CGSC-behavior oc t 2L while \Ii\ oc t L becomes 
much larger. 
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Fig. 15. L = 6 and W = 1: Amplitudes |7)| (lines with sym- 
bols) and |7t| (lines without symbol) as a function of t, ob- 
tained using a disordered sample with increasing values of U 
and $ — 0.05^o- The thick solid lines give the expected lattice 
behaviors. 
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Fig. 16. L = 9 and W = 1: Amplitudes |7;| (line with sym- 
bols) and \It\ (line without symbol) as a function of t, ob- 
tained using a disordered sample with increasing values of U 
and $ — 0.055&0- The thick solid lines give the expected lattice 
behaviors. 



In Figs. ITol and 1161 the amplitudes \Ii\ and \It\ are 
shown as a function of t for W = 1 and various values of U, 
using disordered samples of size L — 6 and 9 respectively. 
The power laws shown by thick lines give the decay of the 
currents t L , t 2L for L = 6 and 9 respectively, as expected 
for lattice regimes (t < t*) with W = 1. The i 3 -behavior 
expected in the BWM-regime requires a smaller disorder 
to be observed. For t > t* , the behaviors are similar to 
those observed for t = 1 and U < U* . In the lattice LWM 
regime (t — > 0), one can see that |Jt| continues to exhibit 



8 Continuum-Lattice crossovers with disorder 

We have shown in section 21 that I oc (K) unless the 
GS begins to occupy states of high momenta k. There- 
fore, the breakdown of the interaction-invariance of I is a 
consequence of a localization-delocalization transition in- 
duced by U or W in momentum space, unless one has a 
level crossing. Without disorder, a jump in the GS total 
momentum K can be induced by level crossings between 
GSs of different K. In the presence of disorder, these level 
crossings become avoided level crossings. We discuss in 
the following sub-section how avoided level crossings limit 
the interaction-invariance of /, outside the lattice regimes 
where the relation / cx (K) is broken. 



8.1 Level crossings below and momentum 
conservation 

When W = 0, K is conserved and one can follow the 
GS of given K as a function of U . When [7 = 0, the 
GS momentum is K = only if the Fermi shell is full. 
For a square lattice, this corresponds to N = 1, 5, 9, . . .. 
For U w < U < U* , the GS is a Wigner molecule with 
harmonic oscillations and a delocalized center of mass. 
This is what we have shown for N — 3 in section. 03 The 
GS energy 

E (K, r.) « E cm (K) + E rel (r s ) + E coul (74) 

depends on K via the motion of the center of mass 
only. In this continuum Wigner regime, the kinetic en- 
ergy of the center of mass is minimum for K = 0. So 
for N — 1,5,9,..., the GS remains in the sub-space of 
K = when U varies and does not exhibit level crossings 
between GSs of different K. For N = 3, this is differ- 
ent, because the GS momentum K ^ at U — and 
must be zero above U w . This yields a GS-level crossing 
at an interaction Ui c which is necessary smaller than U w . 
In the presence of a weak disorder, this crossing becomes 
weakly avoided. Following the true GS as a function of 
U, we expect that its total momentum will be conserved, 
outside the avoided crossing where K exhibits a sharp 
drop from a value K ^ to K = 0. This is illustrated in 
Fig. El where one can see the change AEq of the GS en- 
ergy induced by <P = O.OS^o as a function of U. AE (K) 
are given when W = and L = 6 for K = (0, 0) and 
2n/L(l, 1). One can see the interaction independent be- 
haviors up to a lattice threshold U*(L = 6) w 100, fol- 
lowed by the lattice decay cx U~ 2 . When we include a 
very weak disorder (W/t = 0.01), and follow the true 
GS, AE (W = 0.01) w AE (W = 0,K = 2n/L(l,l)) 
up to the avoided crossing taking place at an interaction 
Uic « 10 where AEq jumps towards the GS-behavior char- 
acterizing W = and K = 2tt/L(0, 0). 
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6 and t = 1: Change AEq of the GS energy 
= 0.005#o as a function of U for different values 



We now study when the relation / oc (K) is broken in a 
disordered lattice of elastic mean free path I and localiza- 
tion length £ without interaction. To obtain the threshold 
values where the lattice-continuum crossover takes place, 
we use the change AEq of the GS energy instead of /, 
and look when the AEq(U) characterizing the three lat- 
tice regimes described in section HO become of the order of 
AE (U = 0). 



8.2 Lattice-Continuum crossover in the ballistic regime 

(L<1) 

When L < I, AEq sa t without interaction and goes as 
t e f f in the BWM-regime. The condition: 



t 



LH 



gives for the threshold U*/t the same equation 



t ~ L ' 



(75) 



(76) 



than Eq. LSI 



8.3 Continuum-Lattice crossover in the diffusive 
regime (I < L < £) 

When the one particle motion is diffusive without inter- 
action, AE (U = 0) oc I(U = 0) oc hek F l/(mLs 2 ) oc 
gt/L 2 (see refs |27U2fi| ). When AEq becomes smaller than 
AEq(U = 0) for values of U > U s t r i P e, the exponent of 
the decay which can be seen in Figs. El and (curves 
for W = O.lt and It) shows that one enters in the lat- 
tice CGSC-regime. Since g oc (t/(LW)) 2 , the crossover 
threshold V stripe is given by the condition AEq(U = 0) « 
gt/L 2 » AE 0: 



stripe ■) 



t 2L L AL 



W 2 L A Ty2L/3-l[/4L/3' 



(77) 



which gives a continuum-lattice crossover taking place in 
the diffusive regime when: 



U f 



stripe 

1 



1/2 



(78) 



Comparing Eq. [7S1 and Eq. 1721 one can see that one gets 
the same relation for the continuum-lattice crossover in a 
diffusive sample (e.g. W s» t for L = 6 — 9) and for the 
BWM-CGSC crossover in the lattice regime of a weakly 
disordered sample (e.g. W w O.Olt for L — 6 — 9). 



8.4 Continuum-Lattice crossover in the localized 
regime (L > £) 

When all the particles are localized without interaction, 
AEq ps iexp— (L/£) where the localization length £ w 
y/t/W for a large ratio W/t (see ref.|25]). When AEq 
decays as a function of U, we can see in Fig. E| that 
this decay corresponds to the lattice LWM-regime. The 
continuum-lattice crossover in strongly disordered samples 
takes place when 



3L-3 



t exp 



which gives an interaction threshold 





(79) 



(80) 



above which one has the LWM-regime instead of a corre- 
lated Anderson insulator (CAf) for a large size L. 

Let us mention that a complete discussion of the var- 
ious possible lattice regimes requires to consider also the 
lattice regimes yielded by large values of W as we have 
considered the lattice regimes yielded by large values of 
U , since the GS can occupy large momenta k either when 
U or W are large. A complete study will require also to 
determine if the CAf regime is in a continuum limit or in a 
lattice limit, as a function of the disorder strength W and 
L. We do not discuss in more details this issue, though the 
data which we show for large values of W/t are certainly 
not in a continuum regime. 



8.5 Sketch of the phase diagram of the different 
lattice regimes obtained for N = 3 

We have sketched in FigED the different lattice regimes 
characterizing = 3 spinless fermions in the plane 
(U/t, W/t). The shaded part of the (U/t, W/t) plane gives 
the continuum regimes which we will study in more details 
in sectional The non shaded part gives the lattice regimes 
which we have studied. Let us summarize the behaviors 
obtained by increasing U for small disorder strengths W. 

When W — 0, one gets between Uf and Uw a contin- 
uum supersolid molecule between the Fermi system and 
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Fig. 18. Different lattice regimes for N = 3 spinless fermions 
in L x L disordered lattices. The shaded part of the (U/t, W/t) 
plane gives the continuum regimes where I does not depend on 
U /t. The thick dashed lines give the threshold values W s tri P e 
and Wioc and W g i a as separating different lattice regimes. U* , 
Uw and Uf are defined in section^ The continuum and lattice 
CAI regimes occuring at large disorders are not separated. 



the continuum Wigner molecule before having a lattice 
Wigner molecule above U* (see section 01 . 

When W/t < 1 and L w 6 — 9, the continuum ballistic 
regimes are followed by three lattice regimes taking place 
above U*/t: a BWM regime for U* < U < U str ipe] a 
CGSC regime for U str ipe < U < Ui oc and eventually a 
LWM regime above Ui oc . The disorder is irrelevant below 
Ustripe , unless the system is at the vicinity of an avoided 
level crossing. 

When W/t is large enough to yield Anderson local- 
ization inside a small scale L w 6 — 9, the Fermi glass 
with Anderson localization (CAI regime) becomes a highly 
correlated solid as U increases, through an intermediate 
regime studied in Ref. [3]. Three examples are given in 
Fig. U3 for W/t = 10,20,40, showing that when U ex- 
ceeds U glass , a localized regime where / is independent of 
U gives rise to another localized regime where I becomes 
independent of W. In those three cases, /; changes its sign 
around U g i ass . As pointed out in Refs. the sign of 

fluctuates from sample to sample for weak U while it 
becomes non-random at large U. 

We study in the following section the effect of U upon 
the continuum limit of a diffusive sample where W = t = 1 
for L — 6 and 9. 



sample fluctuations. This avoided crossing yields a change 
of the sign of //, a sharp drop of |/;| and singularities in 
\It/h\ (see Figs. 1111 and 11211 . Since the total current I is 
independent of U/t outside level crossings in the contin- 
uum limit, one needs to study the local currents and the 
corresponding particle densities for detecting the different 
continuum regimes. 



9.1 Continuum regime for a Coulomb Guided Stripe of 
Current 

As one can see in Figs. ITD and H*2l \It/h\ — 1 below U*/t, 
while /; and It are independent of U instead of exhibiting 
the lattice decays. This shows us that the regime of cur- 
rent stripes is robust, and persists outside the lattice limit, 
in the continuum diffusive limit. This continuum CGSC 
regime is illustrated in Fig. The left figure shows the 
map of local currents obtained in a disordered sample with 
N = 3, L = 9, W = 1, t = 1 and U = 50. The current 
exhibits the Coulomb guided flow along the diagonal di- 
rection, while one can see in Fig. ^| that I does not vary 
as a function of U around U = 50. The corresponding 
site occupation numbers given in Fig. 1191 (right) form an 
extended diagonal stripe, instead of being localized upon 
three main sites, as in the CGSC lattice regime shown in 
Fig. □(right). 
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Fig. 19. Coulomb guided stripe of current in the continuum 
limit. Disordered sample with 7V = 3, L — 9, W = 1, t = 1 
and U = 50. Left figure: map of currents for $ = 0.05^o- Right 
figure: corresponding site occupation numbers n? 



9 Analysis of the continuum regimes of a 
diffusive sample 

For a ratio W/t ps 1, there is a diffusive motion when 
U/t = and L w 6 or 9. We study the continuum limit 
where / does not exhibit the power law decays charac- 
teristic of lattice regimes. For values of U/t smaller than 
the lattice threshold U*/t, one can indeed see in Figs. El 
and EH that \Ii\ and \I t \ are almost independent of U, ex- 
cepted around the value U\ c where an avoided level cross- 
ing occurs. Uic depends on W and t and exhibits sample to 



9.2 Diamagnetic or Paramagnetic currents? 

Before considering weaker values of U to be closer to 
the quantum Fermi limit, it is useful to mention a the- 
orem proved by Leggett |29j for one dimensional spinless 
fermions with arbitrary interaction and external poten- 
tials. This theorem allows to determine the parity of the 
number of particles carrying / from the sign of /. Using 
a certain variational ansatz for the GS wave-function, one 
can show that I is diamagnetic for an odd number N of 
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particles when ss n<?o or for an even number N of par- 
ticles when <P « (n + l/2)^>o> where n is an integer. This 
means that the charge stiffness 



D = 



(-!)"£ 



(£ (0) - Eofr)) > 



(81) 



in one dimension. To odd numbers N correspond diamag- 
netic currents, while they are paramagnetic if N is even. 
This result has always been verified [3U] in numerical stud- 
ies of one dimensional systems with arbitrary disorders 
and interactions. The argument is based on the study of 
the nodes of the wave functions, a simple exercise in one 
dimension which becomes highly nontrivial and still un- 
solved in higher dimensions. If the transverse dimensions 
are small, one can argue that nodal surfaces having a two 
dimensional topology will cost too much kinetic energy, 
and that the sign of / is still given by the one dimen- 
sional rule. In the square lattices which we consider, we 
do not use Leggett's rule because the system is a narrow 
stripe, but because the Coulomb repulsion gives rise to 
a narrow stripe of current in the 2d lattice, reducing the 
two dimensional dynamics to a simpler one dimensional 
contribution. From the sign of /, one can know the parity 
of the effective number of particles which give rise to the 
observed one dimensional stripe of current. 
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Fig. 20. Behavior obtained near the avoided level crossing 
{Uic) for a disordered sample with N = 3, L = 9, W = 1, t = 1 
and U = 15. Left figure: map of currents for $ = 0.05^o- Right 
figure: corresponding site occupation numbers n?. The n- axe 
similar to those obtained in the same sample with U = 50 
(see Fig. I19H while the currents is the superposition of two 
perpendicular diagonal motions. 



9.3 Disordered supersolid regime 

When we continue to decrease U, to eventually get a Fermi 
glass of three independent particles, we observe for N = 3 
a level crossing around t// c where the sign of / changes. 
/ is diamagnetic above Ui c and paramagnetic below Ui c . 
When U is not too weak, the local currents forms a stripe 
of an almost one dimensional topology. This allows us to 
use Leggett's rule to argue that a diamagnetic current 
stripe is carried by N = 3 particles, while the paramag- 
netic current stripe which persists below Ui c is due to a 
single delocalized pair (N — 2 even) in the background of 
an almost localized third particle. Below U\ Cl the corre- 
lated pair gives a dominant paramagnetic current, while 
the third localized particle gives rise to a negligible 2d 
current of random sign. This disordered supersolid regime 
seems to be associated with the 1/VF-behavior of I shown 
in Figs. [Bland [H 

Using the same sample which displays the continuum 
diamagnetic stripe shown in Fig. lliJI for U — 50, one shows 
the map of currents and the occupation numbers n9 in Fig. 
ED and Fig. Ell for U « U lc w 15 and U = 7 respectively. 
The n? are close to those obtained in the same sample for 
U = 50 (see Fig. 119(1 . They form approximately the same 
stripe, though its width becomes broader. However, the 
map of currents are very different. As one can see in Fig. 
PHI for U = 7, the current stripe shown in the left figure 
is perpendicular to the density stripe visible in the right 
figure, while there is a superposition of two perpendicular 
diagonal flows near the avoided level crossing ( Fig. EOJ). 
The width of the paramagnetic stripe of current increases 
to give a two dimensional pattern of local currents when 
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Fig. 21. Behavior obtained below the paramagnetic- 
diamagnetic crossover for a disordered sample with N = 3, 
L = 9, W = 1, t = 1 and U = 7 (r s w rf ) . Left figure: map 
of currents for $ — 0.05^o- Right figure: corresponding site 
occupation numbers n? . The left and right figures are now dis- 
connected, I flowing mainly perpendicularly to the axis where 
the ni are maxima. 



9.4 # . ^o/2 and # /3 harmonics of 

To study more precisely how melts the Wigner molecule 
when U — * in a disordered sample where L = 9, 
t = W = 1, we give in Fig. [52/; as a function of for the 
values of U where one expects the end of the Fermi glass 
(U = 5), a disordered supersolid (U — 15) and the con- 
tinuum CGSC-regime (U = 50) respectively. For large U, 
Leggett's rule gives a diamagnetic molecule of N = 3 par- 
ticles, when the dynamics remains inside a Id stripe. This 
rule is indeed observed for U = 50 where /;(<£) exhibits 



Z. A. Nemeth and J.-L. Pichard: Persistent currents 



17 




0.2 0.4 0.6 0.8 1 

3>/<t> 



Fig. 22. /; as a function of <?/^>o for a disordered sample with 
N = 3, L = 9, W = 1 and t = 1 and three values of U. The 
negative (positive) currents are paramagnetic (diamagnetic) 
for an infinitesimal positive flux. One can see a change of sign 
between U = 15 and U = 50. 



an amplitude as large as for weaker interactions with no 
trace of other characteristic period than <Pq. In contrast, 
the curves obtained for U = 15 and {7 = 5 have reacher 
harmonic structures. 

is an odd function of $ = —/(#)) of period 

which can be expanded as: 

h (#) = V2 ^ /„ sin ( 2™— ) (82) 



where the Fourier components /„ are 

7„ = ^ //(^V^sin^Trn^ d$ (83) 

We give in Figs. EU and EU the <£ , #o/2 and <P /3 
harmonics of as a function of r s for two disordered 

samples of size L — 6 and 9 respectively, with W = t = 1 . 
These harmonics, mainly I\ and J2, give the largest con- 
tributions. One can see that they are negative (paramag- 
netic) up to rf. Around rf , I 3 becomes positive, to reach 
a maximum for r s sa 25 where 1\ w 0. Above r s « 25, the 
odd harmonics l\ and 7 3 are diamagnetic, while the even 
harmonic I2 remain paramagnetic. This is what Leggett's 
rule gives if each /„ was due to a n body motion in one 
dimension. This allows us to give for each sample the char- 
acteristic values where the current topology becomes one 
dimensional and Leggett's rule applies. For U ~ 40 in the 
sample shown in Fig. 1241 the diamagnetic Zj is maximum 
while the two others are negligible, and one has the contin- 
uum stripe (CGS-regime) shown in Fig. 1191 The period <Pq 
should be expected for a continuum disordered stripe, a 
L/3,L/3 translation of each of the 3 particles being equiv- 
alent to have a single particle enclosing <Pq. 
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Fig. 23. Harmonics /„ characterizing the period <l>o/n of the 
longitudinal current -/)($) as a function of r s for a disordered 
sample with N = 3, L = 6, W = 1 and t = 1. 
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Fig. 24. Harmonics J n characterizing the period <Po/n of the 
longitudinal current li(^) as a function of r s for a disordered 
sample with 7V = 3, L = 9, W — 1 and t = 1. 

10 Summary 

Considering only N — 3 spinless fermions in a 2d square 
lattice with random potentials, we have obtained many 
different regimes by increasing the strength U of the elec- 
tronic correlations. Studying how / depends on the system 
parameters, we have distinguished the continuum regimes 
from the lattice regimes. Increasing U in a disordered sam- 
ple of size L, we have studied the effect of electronic cor- 
relations in the ballistic, diffusive and localized regimes 
without interaction. Looking at the map of local currents 
induced by a flux we have identified two regimes rem- 
iniscent of phases discussed in other fields, the stripes in 
the theory of strongly correlated electrons and the super- 
solids in the theory of quantum solids. 

In the regime of the current stripe, the observed diag- 
onal motion is clearly characteristic of the oblique shape 
of the Wigncr molecule with N = 3. For larger N, it is 
likely that one has a regime where the current flows along 
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the axes of the Wigner crystal. If those axes form an an- 
gle with the shortest direction enclosing the flux, the ratio 
It/ h 7^ 0. When the electron crystal is pinned by disorder, 
this suggests that one can have an anisotropic resistivity 
tensor inside a certain range of density. Similar behaviors 
are observed in cuprate oxides with high temperature su- 
perconductivity in a certain range of chemical doping and 
in 2DEGs under high perpendicular magnetic fields. More 
directly related to our problem, the existence of stripe 
phases in 2DEGs is also discussed in Ref. |31| . 

The nature of the intermediate regime between the 
solid and the liquid is a more delicate issue. We suggest 
that one could have a "supersolid" . This concept was in- 
troduced long ago to describe the quantum melting of 3d 
quantum solids at very low temperature when the pres- 
sure is varied. Decades of studies of Helium-3 (fermions) or 
Helium-4 (bosons) atoms have not allowed to reach a firm 
conclusion. Very recent experiments by Kim and Chan in 
solid Helium-4 [32 , 33 , 34 , 35J can be interpreted as an ob- 
servation of an apparent superfluid component, suggest- 
ing that a melting solid could have a supersolid fraction 
at certain intermediate pressures. For electrons in two di- 
mensions, a supersolid phase is a possibility mentioned in 
Refs jHlEl- One can also mention other possibilities. For 
instance a quantum liquid crystal |S] (quantum hexatic 
phase) , in analogy with the thermal melting of an electron 
crystal of very large r s . Clearly, more works are needed to 
understand more precisely the exact nature of the 2DEG 
at intermediate values of r s . Let us conclude by summa- 
rizing three facts characterizing the intermediate values of 
r s : (i) the existence of an unexpected 2d-metallic behav- 
ior given by transport measurements, (ii) the particular 
maps of local persistent currents obtained in exact nu- 
merical studies using a few spinless fermions, obtained in 
this work or in Refs 4,5,6,78,9 , (iii) the results of a re- 
cent fixed node Monte Carlo study \H>\ using N w 50 — 200 
spinless fermions and extrapolated to the limit N — > oo 
which show that the nodal structure of a Slater determi- 
nant of delocalized Bloch waves gives a smaller GS energy 
than the trial wave- functions previously used to describe 
either a pure solid or a pure liquid, for intermediate values 
of r s (30 <r s < 80). 
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provided through the European Community's Human Po- 
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